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We consider a supersymmetric hybrid inflation scenario in which the U(l) i?-symmetry is 
explicitly broken by Planck scale suppressed operators in the superpotential. We provide an 
example with minimal Kahler potential, with the i?-symmetry breaking term relevant during 
inflation being aS 4 , where S denotes the well-known gauge singlet inflaton superfield. The 
inflationary potential takes into account the radiative and supergravity corrections, as well 
as the soft supersymmetry breaking terms. For successful inflation, with the scalar spectral 
index in the currently preferred range, n s fa 0.97 ± 0.010, |a| < 10~ 7 . The tensor to scalar 
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ratio r < 10" 4 , while |dn s /dlnfc| ~ ©(lO" 3 ) - 0(1O- 4 ). 
PACS numbers: 98.80.Cq 



^ ■ 1. INTRODUCTION 

Modern cosmology has seen rapid developments due to experiments such as COBE and WMAP. 
Augmenting their unprecedented successes will be Planck, which may for the first time yield direct 
evidence of inflation. At the same time, there are enormous strides being made in particle physics, 
with the LHC having made perhaps the greatest single discovery in the field in decades, while the 
testing of supersymmetry (SUSY) is highly anticipated. For the first time particle and cosmological 
models can be tested with precision, and the deep connections between these two fields motivate 
us to consider the effects that particle physics considerations have on inflationary models. This 
has led to models such as, among others, SUSY hybrid inflation. 
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The standard version of SUSY hybrid inflation |lrl3( remains one of the most successful and 
, well- motivated inflationary models. It is the most general non-trivial model one can write with a 

gauge singlet field S and supermultiplets <3? and $ that respects U(1)r, such that the latter two 
fields belong to non-trivial representations of the gauge group G. It has a connection to particle 
physics in that within it grand unified theories (GUTs) are naturally incorporated. In this model 
the gauge group G is broken to a subgroup H at the end of inflation, and the energy scale at 
which this occurs is related to local temperature anisotropies in the cosmic microwave background 
radiation [lj]; this scale happens to be close to the GUT scale, indicating that G may be related 
to a GUT. In addition, supergravity (SUGRA) corrections remain under control [4j because this 
model can yield solutions within the WMAP nine-year 2a bounds without trans-Planckian inflaton 
field values. In this model, where only (minimal) SUGRA and radiative corrections are added to 
the global potential, a numerical lower bound on the scalar spectral index, n s 0.985, develops. 
This is somewhat disfavored, although it is within the 2a bounds, n s = 0.971 ±0.010 [B]. However, 
the addition of a linear soft-SUSY breaking term [6| reduces n s significantly 0]; this term is also 
important in explaining the MSSM //-problem [8l4l0l|. It has been shown that small gravity waves 
(tensor-to-scalar ratio r) are produced with minimal Kahler (retaining the lowest-order SUGRA 
term), radiative corrections, and positive soft mass squared and negative linear soft-SUSY breaking 
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terms in the potential. Specifically, we know from [?J that one can obtain r < 10~ 5 in the red-tilted 
region in this case, although r is seven orders of magnitude smaller than the upper limit at the 
central n s value. By contrast, the WMAP nine-year upper limit is r < 0.13. Planck will be able 



to detect r of the order ~ 10 lj|; if gravity waves exist at this order of magnitude, Planck will 



greatly reduce the number of viable inflationary models, or at the least, will rule out a considerable 
amount of parameter space in models in which large gravity waves can be produced. 

One potential drawback of this model is that, since symmetry breaking occurs at the end of 
inflation, topological defects, such as magnetic monopoles and domain walls, can be produced in 
quantities sufficient to contradict experiment. This can be remedied by infla ting along a shifted 
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track such that the symmetry breaking precedes observable inflation [l|, [l2j, [l3|]. However, one 
can solve this problem in the standard case with an appropriate gauge group choice in which no 
topological defects are produced, such as G = SU(5) x U(l) y, referred to as flipped SU(5) [bH ]. 

When doing inflation in the context of flipped SU(5) [19. Il6l|. the i?-charges are determined in 
part by the inflationary sector. By imposing a U(1)r symmetry ( "i?-symmetry" ) , as is standard 
practice, one ensures that proton decay occurs only via the six-dimensional operator (the 5^5^ 
term is disallowed). The importance of this is that we preclude rapid proton decay. However, one 
also prohibits terms such as the quartic couplings 10jl0jl0#10#, which give rise to right-handed 
neutrino masses. (Note that the decay of Majorana right-handed neutrinos can explain the observed 
baryon asymmetry via leptogenesis |17||.) If we reassign i?-charges such that 10jl0jl0# 10h is 
allowed, we end up prohibiting the Yukawa term which gives rise to down- type quark masses 
(yij^ 10;10 j5h)- This outcome is not necessarily catastrophic, since the relevant quark masses 
may be generated radiatively. On the other hand, one may invoke a "double seesaw" mechanism 
[181 ] to account for the lack of large right-handed neutrino masses. A simpler solution to this 
problem is to allow higher-order (Planck scale suppressed) i?-symmetry violating terms in the 
superpotential, while enforcing i?-symmetry for renormalizable terms. With this motivation in 
mind, we wish to explore the effects of these additional terms on the inflationary dynamics. 

In this letter, we consider the inflationary implications of such a scenario, within the context of 
standard SUSY hybrid inflation, which we briefly review in the Section 2. In Section [3] we discuss 
our model including the leading nonrenormalizable terms, and in Section 0] we detail our results. 
There, we compare our results with WMAP nine-year data, and discuss our model in light of the 
Planck mission. 



2. REVIEW OF THE STANDARD HYBRID INFLATION MODEL 

The most general non-trivial renormalizable superpotential one can write involving a singlet 
superfield S and two conjugate supermultiplets $ and $ that preserves a gauge group G and 
U(1)r i?-symmetry is [H, 0] 

W = kS($<S> - M 2 ), (1) 

where M is the energy scale at which G breaks and k is a dimensionless coupling which we take 
to be positive without loss of generality since we can absorb the phase of k into that of S. The 
global SUSY F-terms are given by 

i 

Here, 2, G {s,4>, (/)}, where s, <fi, and 4> are the scalar components of the superfields S, and <l>, 
respectively. We choose to set the D-terms to zero by imposing |<3?| = |<3?|, for convenience. 
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FIG. 1. The tree-level, global scalar potential V in standard hybrid inflation. The variables s and (j> are the 
scalar components of the superhelds S and $. 

Using Equation ([2]), the tree-level global SUSY potential in the D-flat direction is 

V F = K 2 {M 2 -\(p\ 2 f + 2K 2 \s\ 2 \(p\ 2 . (3) 

A plot of this potential in field space is shown in Figure [H Inflation proceeds along the local 
minimum \<f>\ = (the inflationary track), beginning at large |s| (top of Figure [1]). An instability 

occurs at the waterfall point |s c | 2 = M 2 , which is the value of \s\ such that = 
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subscript "c" denotes "critical", and the symbol \s c \ will be used later to denote the dimensionful 
inflaton field at the critical point; we maintain the same notation here for consistency.) At this 
point the field falls naturally into one of two global minima at \(j)\ 2 = M 2 . This coincides with the 
breaking of the gauge group G. At large |s|, the scalar potential is approximately quadratic in \<p\, 
whereas at \s\ =0 Equation ([3]) becomes a Higgs potential. 

Along the inflationary track the potential is flat (V = k 2 M 4 ), and thus one cannot drive 
inflation. Onedoop radiative corrections (RC), which should be added for consistency in any 
case (since SUSY is broken during inflation), can be used to drive inflation. SUSY is restored 
after inflation, when the field evolves to one of its global minima (where V = 0). The radiative 
corrections [19] involve the function 



F(x) = \ 



. a . , (x A — 1) 9 , x 2 + 1 , / K 2 m 2 x 2 
{x 4 + 1) In ^ j— i + 2x 2 In -= + 2 In 



x 4 x 2 -i V Q 



where x = |s|/|s c | is a convenient reparametrization of the inflaton field, N is the dimensionality 
of the representation of the fields $ and and Q is the renormalization scale. 

3. PLANCK SCALE SUPPRESSED i?-SYMMETRY VIOLATION 



We now wish to determine the effects of allowing i?-symmetry violation beyond the renormal- 
izable level. First, let us list which additional terms one can consider in this type of model. The 
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three lowest-order nonrenormalizable ii-violating terms one can write with the aforementioned 
superfields (respecting gauge symmetry) are 

mp mp mp 

where a, /?, and 7 are dimensionless, and are sufficiently small such that each term is a perturbation 
about the standard case. Further, mp = M P /V8tt « 2.4 x 10 18 GeV is the reduced Planck mass. 
Along the inflationary track, only the first term will lift the potential. Therefore, in this letter we 
consider solely the inflationary ramifications of the S 4 term, so that our superpotential is 

W = kS($$ - M 2 ) + — S 4 . (4) 

mp 

Other inflationary tracks may be available via the inclusion the f3 and 7 terms. Inclusion of these 
terms may lead to a form of shifted inflation; however, we do not discuss this here. 

It is important to ask the following question: Can we drive inflation with this S A term, without 
radiative corrections nor any additional terms? Let us compute the global, tree-level scalar potential 
along the inflationary track. We do this via Equation ([2]), which yields the dimensionless potential 
(V = y/mp 4 ) 



V F 



= K 2 m A - 8cos((9 Q + 36 s )\a\ x 3 Ks c 3 m 2 + 16x 6 \a\ 2 s c 6 , (5) 

0|=O 



where 6 a and 9 S are the phases of a and s, respectively, and where we have defined the following 
dimensionless parameters: 



\s\ _ \s c \ _ M 

|,s~ c |' c mp' mp 



The symbol |s c |, as before, denotes the inflaton field at the waterfall point and its dimensionless 
value s c is given by 

- Km 2 + s 2 (k + 4|a|s c ) = 0. (6) 

It can be shown that, using just Equation ([5]), one cannot obtain a red-tilted spectrum while 
simultaneously satisfying the slow-roll conditions (to be defined somewhat later). An analytical 
calculation reveals that, by imposing the condition n s < 1 in the slow-roll approximation, the 
inflaton field at the start of inflation is necessarily trans-Planckian (see Appendix[2]). We therefore 
cannot, in this scenario, achieve a suitable spectral tilt without additional terms. 

In order that our model yield more experimentally favorable results, we include soft [6(] and 
SUGRA corrections (3, 0] to the global plus RC potential. The soft terms are derived in a gravity- 
mediated SUSY-breaking scenario (2^1; including the soft mass squared terms, they are 

am 3 / 2 Km 2 s c x , m 3 / 2 2 s 2 x 2 , bm 3 / 2 \a\ s c 4 x 4 , 

where the last term is a direct consequence of our S 4 term in W, and m 3 / 2 is t ne gravitino mass 
(~ TeV) divided by the Planck scale. We write the effective coefficients of the soft terms as 

a = 2 [2 cos(e s ) - \A\ cos(6 A + 9 S )} , 

b = 2 [\A\ cos(8 A + 9 a + 49 s ) + cos(6 a + 40,)] , ( ' 



where each 6i, i G {^4, a, s}, is the phase of a complex parameter, and A is the trilinear coupling. 
Note that we cannot take a real without loss of generality, since we have already absorbed the 
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phase of k into that of s; therefore, we consider the most general case where a, s, and A are 
complex. 

While 9a and 9 a are components of couplings, 9 S is a dynamical field. For the sake of simplicity, 
we minimize the potential with respect to 9 S so as to define the inflaton field purely as \s\. As a 
result (see Appendix [1]) , we choose the following values of the phases 

9 S = In , 6a = nir , 9 a = pit, 

such that /, n, and p are all odd integers. With these choices, the effective coefficients are 

a = -2(2 + |4|) , b = 2(\A\-l), 

and in conjunction with these phase choices we additionally impose the condition that \A\ < 1, 
or equivalently b < (see Appendix [T|). Henceforth, we drop the bars on A and a with the 
understanding that they represent the moduli of the corresponding complex quantities. 
The SUGRA scalar potential is given by 

Vf;SUGRA = e K/mp2 (K^D Zi WD z *W* - 3m P " 2 \W\ 2 ) , (8) 

with 

Kij = ^-^ , D Zi W = ^+m p - 2 ^W , D Z *W* = {D Z . W)* . 
J ozidz* ozi dzi 

Note that Equation (|8|) reduces to Equation ([2]) in the limit nip — > oo and where minimal Kahler 
is used. (Minimal Kahler is defined as |S| + |<£| + |3>| .) 

We include SUGRA correction terms up to sixth order in the inflaton field |s|, consistent with 
our inclusion of the a 2 |s| 6 global SUSY term; they are: 

12 

— AC 2 m 4 s c 4 x 4 , — m 4 K 2 s c 6 x 6 , — 12k m 2 a s c 5 x 5 . 

Hence, with the addition of the soft SUSY-breaking, SUGRA, and 1-loop radiative correction terms 
to Equation (0), the full scalar potential, scaled by 1/mp 4 , becomes: 



rn 



V = K 2 m A -8an s c 3 m 2 x 3 + 16 a 2 s c 6 x 6 + — " F(x) 

8tt 2 

+ am 3 / 2 m k s c x + bm 3 / 2 a s c x +m 3 / 2 s c x (9) 

1 2 

H — m 4 fi; 2 Sf. 4 x 4 H — m 4 K 2 s r e x e — 12k m 2 a s r 5 x 5 . 

2 3 

In solving the essential cosmological equations, we employ the slow-roll approximation through- 
out, in which inflation occurs while the slow-roll parameters are less than unity. In our notation 
these are written as: 

1 fV\ 2 l V" , 2 _ 1 V"V 



4s c 2 V V / ' ' 2s c 2 V ' s 4s c 2 V 2 ' 

Here, the prime ( ' ) denotes a derivative with respect to x. Inflation ends either when the slow- 
roll parameters become unity, or when the inflaton field reaches the waterfall point at x = 1. 
Observable inflation starts at xq, defined at the pivot scale fco = 0.002 Mpc" 1 , and ends at x e . 
With this, the number of e-foldings becomes, to leading order, 

2 r° v 



N ^2s c 2 j —dx, (10) 



6 



Fundamental 


Range 


Scale 


Derived 


Constraining range 


parameter 




type 


quantity 




K 


[10- 6 ,1] 


log 


n s 


[0.92,1.02] 


m 


[io-Mo- 1 ] 


log 




« 0.971 ± 4cr 


a 


[io- 14 ,io- 8 ] 


log 


A 2 


[2.271,2.583] x 10~ 9 


A 


[io- 10 ,i] 


log 




= 2.427 x 10~ 9 ± 2a 


X 




linear 


r 


< 0.13 


Sc 


(0,1) 


linear 


N 


[50, 60] 



TABLE 1. These are the ranges specified for the fundamental parameters in Equation ©, and constraints 
placed on derived quantities, that we have used in our numerical calculations. Note that x can take on any 
value between the waterfall point and the Planck scale. The experimental bounds on r, A 2 ^, and n s are 
from the WMAP nine-year analysis + eCMB + BAO + H @. The numerical constraints on the quantities 
r and A 2 ^ are the experimental bounds; however, the numerical constraints on n s differ slightly from the 
experimental bounds. This has been done for ease of plotting. 



while the usual definitions hold for 

r«16e , n s «l-6e + 2T7 , -^r « 16er/ - 24e 2 - 2£ 2 . (11) 

dliifc 

The amplitude of the curvature perturbation is given, to leading order, by 

s 2 V 3 



For higher order expressions see [2l| . Note that Equations ([TT]) and (fT2j) are evaluated at the pivot 
scale. 

In our numerical calculations we take rn. 3 / 2 = 1 TeV/mp, Q = 10 15 GeV/mp, and since we are 
implicitly embedding our model in flipped SU(5), we take N = 10 [16]. In addition, we impose 
the ranges in Table [TJ 



4. RESULTS 
A. Overview 



Previous studies have shown that small gravity waves are generated using minimal Kahler and 
a TeV-scale positive soft SUSY-breaking mass squared term (i.e., m 3 /2 2 £ 2 , with m 3 / 2 ~ 10 -16 ) 
, [l6| . Specifically, when the lowest-order SUGRA correction term and a negative linear soft 
term (a = —1) are added to the global SUSY plus TeV-scale positive soft mass squared plus 
RC potential, one finds that r ~ 1CU 12 ' 5 around the WMAP nine-year central value n s = 0.971 
0| . Alternatively, using non- minimal Kahler in shifted inflation with positive TeV-scale soft mass 
squared and a = 1, 0, or —1 [13] (See also 22-24] for further references.), or non-minimal Kahler 



with with same a values, and large, positive soft mass squared terms (~ 10~ 5 ) 25J, one can generate 
r ~ 10 -2 with red spectral tilt. In this paper we find that the solutions follow curves of a similar 
shape to those presented in 0], as can be seen in Figure El By employing minimal Kahler, positive 
TeV-scale soft mass squared terms, and a negative linear and a negative cv-dependent quartic soft 
term, we obtain in this paper r ~ 10~ 8 ' 5 around n s ~ 0.965 for a = 10~ 9 ; in [3], one obtains r 
values four orders of magnitude lower than this, at n s ~ 0.965. As we will describe, this model 
yields even larger gravity waves (~ 10 -4 , see Figure 5(b)) with red spectral tilt, and we expect 
that with non- minimal Kahler this model can yield solutions similar to 25J. While the full set of 
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FIG. 2. The tensor-to-scalar ratio, r, versus the scalar spectral index, n s , is depicted. Three curves for large 
values of a are shown in 2(a) and a larger range of a is taken in 2(b) Here, the number of c- foldings and 
A have been hxed at 50 and 1CP 4 , respectively. With a = 0, we produce solutions closely matching the 
a = — 1 case in [7]. Note that this curve (in 2(b) ) docs not produce false vacua. 



results is outside the reach of current experiments such as Planck, a model in which solutions are 
tending toward larger r solutions is nonetheless preferred. 



B. The Effect of the Parameters on the Model 



Our potential is dependent upon the inflaton field x and the parameters A, a, n, and m. Our new 
parameter a, which parametrizes the amount of ii-symmetry violation beyond the renormalizable 
level, yields qualitatively and quantitatively distinct results from the standard case. 

The negative a-dependent terms in ([9]) create false vacua in some regions of parameter space, 
i.e. the general behavior of the potential changes from that of Figure |3(a)|to that of |3(b)[ (Note 



V(x) 



V(x) 



X 



X 



(a) 



00 



FIG. 3. The qualitative change in behavior caused by the negative a terms. Figure (a) depicts a potential 
that is well-behaved, i.e., the field will roll toward the global minimum, while Figure (b) depicts a potential 
with a false vacuum. Mathematical solutions which produce false vacua are not acceptable inflationary 
scenarios. 



that the innaton field rolls from right to left in these figures.) We cannot produce a successful 
inflationary scenario from Figure |3(b)| as the system will become trapped in the false vacuum. 
Rejecting solutions for which this occurs produces gaps in theparameter space such as those seen 
in Figures [2] and [71 Note that these vacua did not appear in [7j]. 
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FIG. 4. The number of solutions to A|, = 2.427 x 1CT 9 (a), and both A^ = 2.427 x 10~ 9 and N € [50, 60] 
(b). 



Figure [2] depicts the effects of a on r. (See Figure [7] for further results. Note the similarity of 
these curves to those in [3].) The potential in Equation ([9]) differs from that in [?J by two higher- 
order SUGRA correction terms, an a-dependent quartic soft term, and two global a-dependent 
terms (a is also different). The effect of a is to raise r, particularly in red-tilted regions. This is 
primarily a result of the global term proportional to a. The a = 10~ 9 curve is raised by three to 
six orders of magnitude for 0.92 < n s < 0.98, as compared to the a = case in [3]. 

Our model greatly benefits from the fact that we can, as noted, generate larger gravity waves 
than the standard (a = 0) case. However, a cannot be raised arbitrarily, since we require R- 
symmetry violation to be small. We find no need to impose an upper bound, though, because our 
parameter study yields a numerical upper bound a ~ 10 -7 . This can be understood mathematically 
by noting that |V'| increases faster than |V| as a — > 10 -7 (from smaller a). The e- foldings constraint 
(fTTj]) becomes impossible to satisfy at large a, because its integrand is suppressed by a large V, 
and the limits of the integral can only be marginally altered. 

We find that xq can vary over at least two orders of magnitude until a ~ 10 -8 ; then, xq is 
compressed to ~ 10. Likewise, the end of inflation is pushed toward waterfall, x e = 1, as a — > 10 



-7 



Thus the distance in x over which inflation occurs approaches an approximately constant value. 

9 (Figure 4(a) ) and then impose the constraint that 



If we take solutions to 



2.427 x 10~ 

the number of e-foldings be between 50 and 60 (Figure 4(b)), we observe that, for many orders 



of magnitude in a, requiring sufficient inflation decreases by at least an order of magnitude the 
number of solutions relative to those obtained merely from the curvature perturbation constraint. 
Achieving a sufficient amount of inflation severely limits the number of viable solutions generally, 
but is most limiting at large a. The curvature perturbation constraint also has a limiting effect as 
a increases, although notice that without ensuring a sufficient amount of observable inflation, one 
can obtain "solutions" up to a = 10~ 5 (Figure 4(a)). 

Numerically, we obtained our results using two independent methods: a continuation method 
and a parameter study. The results of the former are seen in Figures [2] and [71 The latter results 



are presented in Figures [5] and [U Figure 5(a) shows the presence of qualitatively new "horizontal 
solutions" in the r 



n. 



plane. If we drop all the a-dependent terms in the scalar potential 



m x 



the number of horizontal solutions 



except the global term linear in a (which is —8a k s c 
in this region increases. On the other hand, dropping this term and keeping the others does not 
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n s n s 

(a) (b) 

FIG. 5. (a) The tensor-to-scalar ratio, r, versus the scalar spectral index, n s , for a — 10~ 9 , A — 10~ 4 , and 
Nq g [50,60]. One can see the "horizontal solutions", which yield r values somewhat above 10~ 8 . (b) The 
tensor-to-scalar ratio, r, versus the scalar spectral index, n s for solutions corresponding to the ranges in Table 
[TJ Solutions are color-coded as follows: blue - (10~ 14 < a < 10 -12 ), green - (10 -12 < a < 10 -11 ), yellow - 
(KT 11 < a < 10~ 10 ), red - (lO^ 10 < a < 10~ 9 ), cyan - (10~ 9 < a < 10^ 8 ), magenta - (10" 8 < a < 10~ 7 ). 




a 



(a) (b) 

FIG. 6. (a) The tensor-to-scalar ratio, r, versus a using the constraints in Table [TJ (b) The scalar spectral 
index, n s , versus a using the constraints in Table[TJ Note that in the red-tilted region where a — ¥ 0, solutions 
are sparse. 



produce a viable inflationary scenario. Since the global a term is the most dominant a-dependent 
term, it is primarily responsible for producing these horizontal solutions. 

We find that limits on r naturally arise in our model (see Figure 6(a)). Solutions producing 
r ~ 10~ 4 can be produced throughout the range of a that we have taken, although all of these 
except the solutions near the upper limit of a correspond to a blue-tilted spectrum. Figure |5(b)| 
depicts this behavior; note that the only large-r solutions corresponding to red spectral tilt are 
colored magenta and cyan (meaning that 10" 9 < r < 1(T 7 ). Figure |6(a)| also indicates that 
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the smallest-r solutions that can be produced are increasingly larger as a increases, so that at 
the upper limit of a only r ~ 10~ 4 can be produced. This effect can be understood by using 
the following approximation of the energy scale of inflation: Vo 1//4 ~ n 3 / 2 m 2 / V \ 2 xo 2 a, where 
A = [(3456 A^l^vr 2 )] 1 / 4 . As mentioned, the integrand in Equation (fTOj) is suppressed by the fact 

that |V'| increases faster than |V| as a increases. To compensate for this effect the K 2 m 4 term in V 
increases, raising the numerator in Vo 1 ^ 4 - This prohibits small-r solutions in large-a regions (recall 
that r oc V 1/4 ). 

The upper limit on r of ~ 10 -4 , approximately constant over many orders of magnitude of a, is 

x m 2 (xo — x e ) 2 . The largest m values obtainable 
x e , and therefore limit r below 10 -4 . Solutions 



a consequence of the Lyth Bound [26] r < 0(1O -2 ) 
(m ~ 10 _1 ) correspond to very small values of xo - 
with m ~ 10~ 2 correspond to x$ ~ 10 and x e ~ 
values obtainable in our model. 
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FIG. 7. Here we plot our numerical results in which the number of e-foldings and A have been kept fixed at 
50 and 10 -4 , respectively. The absolute value of the running of n s , log 10 | j ^ * k \ , is plotted in |7(d)[ 



The parameter A, arising from the gravity-mediated soft-SUSY breaking terms, does not have 
any discernible effect on our results in the red-tilted region. This is expected from the fact that 
the soft terms are suppressed by the gravitino mass (~ TeV). 

The plots in Figure [7] depict the effects of increasing a. Figures 7(a) , 7(c) and |7(d)| are in 
direct reference to [7f. We can see from Figure 7(c) | that large a boosts k, especially in the red- 



tilted region. Similarly, Figure 7(a) shows that, for the same breaking scale M, k is boosted as a 



increases. Note that m is also boosted by a (Figure [7(b) [ ). We have discussed that larger gravity 
waves are produced via increasing a. This effect is evident from Figures [7(b)1 and [7(c)} which 
show that m and k values are raised as a increases, particularly in red-tilted regions. Figure [7(d) | 
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depicts the absolute value of the running of n s , log 10 | d d t " s fc | . In our model, is negative, and 



I dn a I 
I dlnfe I 



O(10" 3 ) -0(10" 



5. CONCLUSION 

In this letter we have considered the effects of the i?-symmetry violating -^S A term in the 
standard hybrid inflationary scenario. By allowing for i?-symmetry violation beyond the renormal- 
izable level, we can give masses to right-handed neutrinos and down-type quarks in a simple way 
within flipped SU(5), providing a particularly well- motivated inflationary scenario. We compare 
our results to the WMAP nine-year parameters n s and r, finding that we can easily generate red 
spectral tilt values within the observational bounds. This type of model can generate larger gravity 
waves than the standard (a = 0) case; raising a from 10 -14 to 10 -9 corresponds to an increase in 
r from 10~ 12 to 10~ 8 at the WMAP nine- year central n s value 0]. This result has been achieved 
using a minimal Kahler potential, positive TeV-scale soft mass squared terms, a negative linear 
and a negative quartic soft term, and SUGRA correction terms up to sixth order in x. Additional 
interesting results were observed in the parameter study, which yielded a region of qualitatively 



new "horizontal solutions" in the red-tilted region, above the main r — n s curve in Figure 5(a) 
The parameter study yields r values as large as ~ 10 -4 in the red-tilted region. Finally, we observe 
that | a | develops a numerical upper bound of ~ 1CJ- 7 , while |^| < 0(1(T 3 ) - 0(HU 4 ). 
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Appendix 1: Minimization of the Potential With Respect to 6s 

The phase-dependent terms in our potential are, including just the global SUSY and soft SUSY- 
breaking terms: 

V(6) = -ci cos(<5) + c 2 (2 cos(6> s ) - \ A\ cos(V>)) + c 3 (\A\ cos(<5 + V) + cos(6» s + 5)) , 
where the phases and coefficients (the latter are dimensionful) are 

5 = e a + 36 s , ip = e s + e A 

8|a||S| 3 M K 2 2m 3/2 |a||S| 4 

ci = , c 2 = 2m 3/2 M^K\s\ , c 3 = '- . 

mp mp 

We seek to minimize this function with respect to to 6 S , thus we impose the conditions 

dVfl 

— p = 3ci sm(5) + c 2 -2sin(0 s ) + \A\ sin(^ - 4c 3 [\A\ sin(J + ^) + sm{9 s + 5)} = 0, 
dVfl 

j = 9ci cos(5) + c 2 [-2 cos(6» s ) + \ A\ cos(V?)] - 16c 3 [|A| cos(5 + V) + cos(6 s + 5)] > 0. 

d6/ 
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The first condition is met trivially by setting 



sin(0 s ) = 



sin(^) = 



sin(5 + tp) = 



sin(0 s + 5) = 0, 



which implies that each argument is an integer multiple of ir. Finally, the second condition can be 
satisfied with 

cos(6 s ) < , cos(^) > , cos(<5) > , cos(6 s + 5) < , cos(5 + ip) > 0. 

The last inequality is a consequence of the previous four, but we can nonetheless ensure that the 
second condition is satisfied by declaring \A\ < 1, which, from Equation (J7J), is equivalent to b < 0. 
With these phases the linear coefficient is a = —2(2 + < 0. Finally, the only phase-dependent 
SUGRA correction term we are including is VgjsuGRA = — ^"^jjf • cos (5), which will always be 
negative with our chosen phases. 



It is interesting to investigate whether inflation can be driven exclusively using the new S 
contribution without extending to SUGRA. In global SUSY, the scalar potential acquires two new 
terms (see Equation ([5])) which impart a nonzero slope along the inflationary valley. For \a\ small 
enough to constitute a perturbation, the |s| 6 term is strongly subdominant to the |s| 3 term, and 
thus the latter is chiefly responsible for meaningful alterations to the inflationary dynamics. 

Some care must be taken in choosing the phases of the complex quantities in this version of 
the model. The new term in the superpotential induces an additional SUSY vacuum appearing 
along the |0| = direction. This vacuum is gauge symmetric and must be avoided if the symmetry 
breaking structure of the model is to be preserved. If the phase factor cos(# a + 36 s) is positive, a 
stable inflationary valley along \(f>\ = leads to this vacuum. Therefore we require this phase factor 
to be negative. To achieve this, one (but not both) of 9 a ,6s must be an odd integer multiple of 
7T. For simplicity, we choose 0g = 0, a = nn (n odd). Additionally, these phases result in a slight 
(~cubic) uplifting of the potential, consistent with a perturbation to the original hybrid model. 

With these phases chosen, we examine the constraint imposed on the parameters of the model 
by the slow roll conditions eo < 1, rjo < 1 . We find that the constraint from 770 is convincingly more 
stringent than that from eo, and so we will concentrate on the former. In terms of dimensionful 
quantities, we may write analytically 



The denominator of this expression is strikingly close to Equation ([6]) specifying the waterfall point 
s~ c , after adjusting for phase differences (which amounts to a sign flip for the a-dependent term). 
Approximating |sq| ~ s c , we have 



Appendix 2: Blue-Tilt in Standard Global Plus S 4 Case 



7/0 = 247np|a||so 



kM 2 + W\a\\s \ 3 /mp 



(kM 2 +4\a\\s \ 3 /m P ) 2 



2Amp\a\ \sq 



kM 2 + 4|a||s | 3 A"-p' 



24mp|a| |s 



(21) 



770 ^ 




Taking 770 < 1 and rearranging, we may cast the slow roll condition in the form 
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Next, we investigate the condition for obtaining a red-tilted spectrum in this model. Since 
t]q > in the present case, n s < 1 requires 6eo > 2tjq according to Equation (fTTj) . or 

-4 

Vo 3 

Using the same approximation stated above, we may write 



■so 



eo - 144 - 

. K / \ S 



Combining this with Equation ([22j) yields a red tilt condition 



\a\ I |£o| 



rap J 18 



— > io- 1 



m P \ f Jc_ 

s ~c J Vl s o| 



(23) 



Taking 1/24 ~ 10 1 in Equation (j22j) . the slow roll and red tilt conditions may be consolidated 
into one expression, 



10" 



V s c 7 Vl s o 



\mp 



so 



(24) 



The middle portion \a\/n may safely be omitted, after which several common factors may be 
canceled from the inequality. Simplifying Equation (124j) readily leads to 

(25) 



I s | > rap- 



In other words, the global-SUSY version of the model under the slow roll approximation may 
only lead to n s < 1 for trans-Planckian values of the inflaton field. Since the SUGRA sector is 
suppressed by powers of \sq \ /mp, it is then inconsistent to treat the model only using global SUSY. 

In addition to these arguments, the assumptions placed on the phases for this version of the 
model are also problematic. Above, we commented that the phase factor cos(9 a + 39s) must 
be negative in order to preserve gauge symmetry breaking. Referring to Appendix HJ however, 
the potential is minimized with respect to 9s for cos(9 a + 39s) > 0. Therefore an unattractive 
choice must be made: allow the phase #s to vary dynamically (potentially ruining inflation in the 
process) , or sacrifice gauge symmetry breaking (a defining characteristic of hybrid inflation) . This 
issue compounds the difficulties in obtaining a red-tilted spectrum via slow roll. Thus we conclude 
that successful inflation may not be implemented in our S 4 model using only global SUSY. 
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